May 10, 2001; accepted by The Astronomical Journal. 
Preprint typeset using lAT^^ style cmulateapj v. 26/01/00 



O 
O 
(N 



(N 
(N 



> 
O 

m 

o 



o 



CENTRAL STRUCTURAL PARAMETERS OF EARLY-TYPE GALAXIES AS VIEWED WITH 

SwARA Ravindranath^'^, Luis C. Ho^, Chien Y. Peng'', Alexei V. Filippenko^, 
AND Wallace L. W. Sargent^ 
May 10, 2001; accepted by The Astronomical Journal. 

ABSTRACT 

We present surface photometry for the central regions of a sample of 33 early-type (E, SO, and SO/a) 
galaxies observed at 1.6 /im (iJ band) using the Hubble Space Telescope (HST). Dust absorption has 
less of an impact on the galaxy morphologies in the near-infrared than found in previous work based 
on observations at optical wavelengths. When present, dust seems to be most commonly associated 
with optical line emission. We employ a new technique of two-dimensional fitting to extract quantitative 
parameters for the bulge light distribution and nuclear point sources, taking into consideration the effects 
of the point-spread function. Parameterizing the bulge profile with a "Nuker" law (Lauer et al. 1995), 
we confirm that the central surface-brightness distributions largely fall into two categories, each of which 
correlates with the global properties of the galaxies. "Core" galaxies tend to be luminous ellipticals 
with boxy or pure elliptical isophotes, whereas "power-law" galaxies are preferentially lower luminosity 
systems with disky isophotes. The infrared surface-brightness profiles are very similar to the optical, 
with notable exceptions being very dusty objects. Similar to the study of Faber et al. (1997) based on 
optical data, we find that galaxy cores obey a set of fundamental-plane relations wherein more luminous 
galaxies with higher central stellar velocity dispersions generally possess larger cores with lower surface 
brightnesses. Unlike most previous studies, however, we do not find a clear gap in the distribution of 
inner cusp slopes; several objects have inner cusp slopes (0.3 < 7 < 0.5) which straddle the regimes 
conventionally defined for core and power-law type galaxies. The nature of these intermediate objects 
is unclear. We draw attention to two objects in the sample which appear to be promising cases of 
galaxies with isothermal cores that are not the brightest members of a cluster. Unresolved nuclear point 
sources are found in ~50% of the sample galaxies, roughly independent of profile type, with magnitudes 
in the range m^'' = 12.8 to 17.4 mag, which correspond to Af™'' = -12.8 to -18.4 mag. Although 
the detection rate of compact nuclei seems favored toward galaxies spectroscopically classified as weak 
active galactic nuclei, wc find no significant correlation between the near-infrared nuclear luminosities 
and either the optical emission-line luminosities or the inferred black-hole masses. 

Subject headings: galaxies: active — galaxies: elliptical and lenticular, cD — galaxies: nuclei — 
galaxies: photometry — galaxies: Seyfert — galaxy: structure 



I. INTRODUCTION 

Prior to the advent of the Hubble Space Telescope 
(HST) ground-based studies of luminous elliptical galax- 
ies showed that the surface-brightness profiles had a cen- 
tral core in most cases, but the effect of atmospheric see- 
ing (typically ^1") made it difficult to discriminate be- 
tween truly resolved isothermal cores and unresolved cores 
(Schweizer 1981; Kormendy 1985a). Kormendy (1985a, 
1985b) studied elliptical galaxies and spiral bulges using 
images taken under excellent seeing conditions (0'.'2-0'.'5) 
and confirmed that isothermal cores were indeed very rare. 
However, some ellipticals showed evidence for isothermal 
cores, and these were generally the brightest galaxies in 
rich clusters. Attempts were made to relate the core prop- 
erties with the global properties by accounting for the ef- 
fects of seeing through image deconvolution (Lauer 1985a). 
However, resolution approaching O'.'l would be required to 
resolve cores in low-luminosity galaxies for which the cor- 



relation between luminosity and core radius implies small 
core sizes (Kormendy 1985b). 

Results from l/-band imaging of early-type galaxies us- 
ing the Wide Field Planetary Camera (WFPC) showed 
that traditional functions used to fit ground-based surface- 
brightness distributions, such as King (1966) or de Vau- 
couleurs (1948) r^/^ profiles, do not provide adequate fits 
for the central (r ^1") regions (Crane et al. 1993; Fer- 
rarese et al. 1994; Forbes, Franx, & Illingworth 1995; 
Lauer et al. 1995). The King models used for giant ellipti- 
cals have central cores with constant luminosity densities, 
which cause the brightness profiles to appear flat in the 
center. However, the HST sluAies argued against the exis- 
tence of such isothermal cores based on the non-zero cusp 
slopes seen even within r <^1". It was evident that some 
galaxies have profiles that can be represented by a single 
power law all the way to the HST resolution limit, while 
others require double power laws, with the inner slope in- 
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terior to some characteristic radius being much shallower 
than the outer slope. This led to the formulation of an 
empirical function, essentially a double power law, to de- 
scribe the surface-brightness profiles. The "Nuker" law 
(Lauer et al. 1995; Byun et al. 1996) has the functional 
form 



/(r) = 2('5-T)/«/fc 



n 



1 



n 



(1) 



where /3 is the asymptotic slope as r ^ oo, 7 is the asymp- 
totic slope as r 0. ri, is the break radius at which the 
outer slope /3 changes to the inner slope 7, a controls the 
sharpness of the transition between the inner and outer 
slopes, and 7^ is the surface brightness at rt- 

Two classes of early- type galaxies can be identified based 
on the value of 7 (Lauer et al. 1995). "Core" galaxies 
have surface-brightness profiles which exhibit a significant 
flattening within a well resolved rt (^0'.'2), such that 7 
^0.3, whereas the profiles of "power-law" galaxies do not 
show any significant break but continue to be steep (7 
^0.5) up to the resolution limit of r ~ Qf.'l. There was 
an apparent dichotomy in the distribution of 7 values. 
Core galaxies have < 7 < 0.3, and power-law galax- 
ies show 7 > 0.5; none seemed to have 7 in the range 
0.3 to 0.5 (Faber et al. 1997). Core profiles are found 
predominantly in luminous, slowly rotating systems with 
boxy or pure elliptical isophotes, while power-law profiles 
occupy less luminous, rapidly rotating galaxies with disky 
isophotes (JafFc et al. 1994; Byun et al. 1996; Faber et al. 
1997). Carollo et al. (1997) obtained Wide Field Plane- 
tary Camera 2 (WFPC2) F-band and 7-band images for 
a sample of elliptical galaxies with kinematically distinct 
cores. Parameterizing the surface-brightness profile using 
the Nuker law, they too found that fast-rotating, disky 
galaxies have steep inner slopes while slow-rotating, boxy 
systems tend to have shallow inner slopes. In spite of these 
trends with global galaxy properties, which are similar to 
those of kinematically normal ellipticals, they suggest that 
the inner slopes may vary continuously between core and 
power-law galaxies. Recently, a detailed analysis of the 
surface-brightness profiles for early-type galaxies has been 
carried out by Rest et al. (2001) using i?-band WFPC2 
images. Their results, also based on Nuker fits to the pro- 
files, present further evidence for a continuous distribution 
of inner slopes. 

This paper analyzes i7-band (1.6 /xm) images of a sam- 
ple of 33 early- type (E, SO and SO/a) galaxies observed 
with NICMOS on HST. Our primary aim is to param- 
eterize the intrinsic distribution of the bulge light. We 
achieve this using a two-dimensional (2-D) decomposition 
technique which deblends the bulge from a central nu- 
cleus and properly accounts for the NICMOS point-spread 
function (PSF). As is well known from the WFPC and 
WFPC2 studies, even early-type galaxies contain consid- 
erable amounts of dust in their centers which corrupts 
the surface-brightness profiles and isophotal parameters 
(van Dokkum & Franx 1995; Verdoes Kleijn et al. 1999; 
Tomitactal. 2000; Tran et al. 2001). Near-infrared (NIR) 
images are better suited for isophotal analysis because of 
the reduced sensitivity to dust extinction and emission 
from young star clusters. They also trace more faithfully 
the underlying stellar population dominating the mass. 



The paper is organized as follows. Section 2 describes 
the sample and data reduction. Section 3 summarizes our 
methods of isophotal analysis, including an introduction of 
our 2-D modeling technique. The principal results concern 
small-scale structures in the inner regions (§ 4), surface- 
brightness profiles and central-parameter relations (§ 5), 
and unresolved nuclei (§6). We discuss the implications 
of our findings in § 7 and present a summary in § 8. Ap- 
pendix A gives noteworthy details on individual objects. 
Throughout this paper distance-dependent quantities are 
based on a Hubble constant of Hq = 75 km s~^ Mpc~^. 

2. THE SAMPLE AND DATA REDUCTION 

Our sample is based on galaxies selected from the Palo- 
mar study of nearby galaxies, a ground-based optical spec- 
troscopic survey of a nearly complete sample of 486 galax- 
ies with Bt < 12.5 mag and declinations greater than 
0° (Filippenko & Sargent 1985; Ho, Filippenko, & Sar- 
gent 1995, 1997a, 1997b). Our study focuses on the 
subset of early-type galaxies with NICMOS images avail- 
able in the HST data archive. This comprises 33 objects 
(14 E, 16 SO, and three SO/a) from different snapshot 
programs. The global properties of the sample galax- 
ies are given in Table 1, along with their nuclear spec- 
tral classifications, which includes five Seyferts, 15 low- 
ionization nuclear emission-line regions (LINERs), three 
"transition" (LINER/H II) objects, one H II nucleus, and 
nine absorption-line nuclei (see Ho et al. 1997a for a de- 
scription of the classification system) . 

The data used in this work consist of images obtained 
in the F160W filter {H band) using the NIC2 and NIC3 
cameras. The NIC2 images have a field of view of 19'.'2 x 
19'.'2 and a pixel scale of 0'.'076; the gain and read noise are 
5.4 e~ count"^ and 30 e~, respectively. The NIC3 camera 
provides a field of view of 51'.'2 x 51'.'2 and a pixel scale 
of 0'.'2; its gain and read noise are 6.5 e~ count"^ and 30 
e", respectively. The majority of the data were taken as 
part of snapshot survey programs, and therefore the expo- 
sure times were relatively short, ranging from 15 to 320 s, 
typically ^150 s. 

Data reduction begins with the standard ("calnicA") 
pipeline processing performed at the Space Telescope Sci- 
ence Institute (STScI; Bushouse et al. 1998). The cal- 
nicA task works on the raw science data and removes the 
instrumental signatures through bias subtraction, dark- 
current correction, nonlinearity correction, and flat fleld- 
ing. A ramp-fitting procedure removes most of the cosmic 
rays from the images. A few additional processing steps, 
as recommended in the NICMOS Handbook (Dickinson 
1999), are required to remove residual image anomalies. 
The NICMOS arrays have a number of bad pixels, some 
of which have very low response and appear dark (cold pix- 
els) and others which have very high or erratic dark cur- 
rent and appear bright (hot pixels) in the images. In addi- 
tion to bad pixels, there are regions of reduced sensitivity, 
termed "grots," which result from flecks of anti-reflective 
paint that have scraped off from the optical baffles. The 
best way to eliminate the effects of bad pixels would be 
to use dithered images. Since the archival images we use 
originate from various observing programs, dithered im- 
ages are lacking for some galaxies, and in a few cases only 
single exposures are available. Therefore, we created a 
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mask of all bad pixels from the flat-field images using the 
ccdmask task in IRAF^ and them used the fixpix task to 
interpolate over the masked pixels. 

The NIC2 camera has a coronographic hole which is 
used when observing targets close to bright objects. In 
the non- coronographic F160W observations used for this 
study, the hole appears in the images as a bright circu- 
lar patch with positive intensity due to excess background 
emission from warmer structures located behind it. The 
position of the coronographic hole with respect to the de- 
tector coordinates is known to have moved with time, and 
this causes two patches to appear in the calibrated image, 
the second one arising from the use of reference files taken 
at a different time. We masked the coronographic hole 
during further analysis. 

The calnicA-processed images contain a residual offset 
or "pedestal" which appears as a time-variable bias that 
varies from one quadrant to another. This variable quad- 
rant bias is believed to result from subtle temperature 
changes in the electronics or the detectors themselves. The 
four quadrants of each NICMOS detector have separate 
amplifiers, thereby resulting in a variation of the pedestal 
among them. The effect of the residual bias, which ranges 
from 0.1 to 0.35 count s^^, becomes evident when the dark- 
subtracted images are multiplied by the flat-field reference 
files, whereby the relative pixel sensitivities dominate the 
images (Boker et al. 1999). We used software developed by 
R. van der Marel at STScI to remove the pedestal effect^. 

An important issue when seeking information at the 
highest spatial resolution is to understand the properties 
of the PSF. The PSF often defines the resolution and sensi- 
tivity limits of the observation and can change with time, 
wavelength, position, and the camera used. At 1.6 fiia, 
the PSF is critically sampled for the NIC2 camera (core 
FWHM « 2.3 pixels, or 0'.'17), while it is undersampled for 
the NIC3 camera (core FWHM « 1.1 pixels, or 0'.'22). In 
NIC2 images the drift in the cold masks with time causes 
variations in the spider patterns and diffraction rings of 
the PSF. However, the PSFs of all the NICMOS cameras 
can be well modeled by the Tiny Tim software (Krist & 
Hook 1999), which takes these effects into account. In the 
absence of PSFs empirically derived from observations of 
bright stars, synthetic PSFs generated by Tiny Tim serve 
as good alternatives for use in photometry, deconvolution, 
and image modehng (Krist & Hook 1997). 

3. SURFACE PHOTOMETRY 
3.1. Isophotal Analysis 

Galaxies appear relatively smooth in NICMOS images 

compared to the optical images, and thus are well suited 
for deriving surface-brightness profiles. We used the el- 
lipse task in STSDAS to perform surface photometry. This 
task assumes that the isophotes of a galaxy can be repre- 
sented by ellipses. The ellipse fit is performed by providing 
as input parameters an initial guess value for the galaxy 
centroid, the ellipticity, and the position angle, and then 
allowing all the parameters to vary with increasing semi- 
major ax;is radius. The initial coordinates for the galaxy 
center are estimated using the imexamine task. However, 



the isophotes of a galaxy are often not perfect ellipses. 
The deviation from the fitted ellipses is quantified by the 
higher-order coefficients in the Fourier series 



N 

E 

n=l 



[A„sin(n^) -|- i?„cos(n^)] , 



(2) 



where N is the highest order fitted and is the angle mea- 
sured counter-clockwise from the major axis of the ellipse. 
The amplitudes of the second-order Fourier terms {A2 and 
B2) are used to obtain the position angle (P. A.) and ellip- 
ticity (e) of the best-fitting ellipse. A perfectly elliptical 
isophote is completely described by the first and second- 
order Fourier terms of the above equation. Non-zero coef- 
ficients for the higher-order terms in the expansion series 
carry information on the shape of the isophotes. These 
intensity coefficients are divided by the semi- major axis 
length and local intensity gradient to measure actual radial 
deviations from perfect ellipses (Jedrzejewski 1987). The 
third-order terms have been found to have significant non- 
zero values when the isophotes get skewed by the presence 
of dust (Peletier et al. 1990). The most dominant mode 
which carries information about the isophote shapes is the 
fourth-order cosine term, and its amplitude B4 is positive 
for disky isophotes and negative for boxy isophotes. A 
non-zero value for the corresponding sine term indicates 
rotation or misalignment relative to the major axis of the 
ellipse which could result from the presence of dust or pro- 
jection effects (Franx, Illingworth, & Heckman 1989). 

The results of the isophotal analysis are shown in Fig- 
ure 1. Note that our ellipse fits are performed on the 
observed (without deconvolution) images. (This is in con- 
trast to the 2-D fits in § 3.2, where we properly account for 
the PSF.) The effects of the PSF can be seen in the lower 
three panels of Figure 1, where interior to r ^0'.'2-0'.'3 
the values of e, P. A., and B4 have large formal error bars 
and often undergo erratic variations. The large fiuctua- 
tions at small radii are also partly caused by the discrete 
sampling and sub-pixel interpolation implemented in the 
ellipse routine (see discussion in Rest et al. 2001). We 
consider these data points to be unreliable. The radial 
variations in ellipticity and P.A. at large radii are caused 
by the presence of features such as dust lanes, weak disks, 
or nuclear bars, and thus reveal morphological peculiari- 
ties in galaxies. In some cases, the variations in ellipticity 
and P.A. may also reflect triaxiality in galaxies. The H- 
band magnitudes were computed using the photometric 
keywords provided in the image header for converting ob- 
served counts to the Vega magnitude system (or the John- 
son if-band magnitude). The conversion is uih = -2.5 
log(count s~^) -|- rrizpt, where m^pt = 21.75 for NIC2 and 
21.50 for NIC3. 

3.2. 2-D Modeling 

Ground-based studies of the surface-brightness profiles 
of galaxies traditionally employ 1-D analysis, whereby the 
radial profile is first derived through isophote fitting, fol- 
lowed by decomposition with a combination of an exponen- 
tial disk and a de Vaucouleurs bulge profile (e.g.. Freeman 
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1970; Kormendy 1977; Boroson 1981; Kent 1985; Baggett, 
Baggett, & Anderson 1998). However, in recent years the 
importance of using the entire 2-D image to obtain more 
accurate surface photometric parameters has been exten- 
sively discussed in the literature (Byun & Freeman 1995; 
de Jong 1996; Moriondo, Giovanardi, & Hunt 1998; Scorza 
et al. 1998; Wadadekar, Braxton, & Kembhavi 1999). 2- 
D image decomposition proves superior to decomposition 
in 1-D when the shape parameters are significant, for in- 
stance when a strong disk is present. 

When decomposing galaxies into constituent compo- 
nents, modeling degeneracy is a serious issue, particularly 
in the context of bulge-to-disk decomposition (e.g., Kent 
1985; Byun & Freeman 1995; Moriondo et al. 1998). For 
example, a combination of the de Vaueouleurs law and an 
exponential disk has been widely used, but they usually fit 
well only over a limited range in radius, and often poorly 
near the nucleus. A certain subjectivity is sometimes in- 
volved, therefore, in deciding where and how to fit — a 
decision which affects the derived bulge and disk parame- 
ters. An additional complication is that there is no reason 
why the popular r^/^ function should be preferred over, 
say, the more general formulation of Scrsic (1968), which 
often may fit just as well, if not better. These problems 
can be serious in 1-D modeling. But 2-D modeling, which 
uses the full spatial information, offers the potential to 
break the degeneracies: bulges and disks can often be de- 
coupled in 2-D because they may have different position 
angles and ellipticities. The appropriateness of the chosen 
functions is also usually more apparent. 

In a similar vein, significant ambiguity can arise when 
attempting to extract compact sources (nonstellar emis- 
sion from active galactic nuclei or nuclear star clusters) 
in the centers of galaxies. The degeneracy here is caused 
by the seeing, which reduces the contrast between the nu- 
clear source and the underlying galaxy profile, which itself 
can be sharply peaked. Earlier studies of nearby nuclei 
extracted central sources by fitting 1-D galaxy profiles si- 
multaneously with a point source (e.g., CaroUo, Stiavelli, 
& Mack 1998). The results can depend sensitively on de- 
tails such as the manner in which the 1-D profile is ob- 
tained. Ferrarese et al. (1994) find that while there is 
no noticeable difference in the surface-brightness profiles 
along the major and minor axes for core galaxies, the same 
is not true for the power-law galaxies. Contamination by 
dust can affect the centering and distort the overall nu- 
clear profile, especially in undersampled images. But in 
2-D, one can again use the shape information present in 
the entire image to one's advantage. Wadadekar et al. 
(1999) illustrate the benefit of using 2-D modeling over 
1-D for extracting central point sources. 

In the context of HST images, Lauer et al. (1998) 
showed that even with WFPC2 images, PSF smearing 
significantly affects a galaxy's profile within (y.'5. To re- 
move the effect one can cither deconvolve the image or 
convolve a model profile with the PSF to match the ob- 
served images. Which strategy to adopt depends on the 
signal-to-noise ratio of the data, and to a lesser extent on 
knowledge of the PSF. Our NICMOS images come mainly 
from various snapshot survey programs and hence do not 
have uniformly high signal-to-noise ratios to enable reli- 
able deconvolution. We also do not have access to near- 



contemporaneous observations of stars to derive empirical 
PSFs. The synthetic Tiny Tim PSFs are not known to 
high accuracy because of thermal stresses and "breath- 
ing" effects. We decided, therefore, that convolution is 
the most transparent way to proceed. 

We parameterize the bulge light distribution of our sam- 
ple of early-type galaxies using a 2-D Nuker profile. This 
was done using a least-squares fitting program, GALFIT 
(C. Y. Peng et al., in preparation), which allows fitting 
a superposition of analytic models (e.g., Sersic/de Vau- 
eouleurs, Nuker, exponential, Gaussian, Moffat) to create 
arbitrarily complex galaxy shapes. The radial shapes of 
the models are generalized ellipses (Athanassoula et al. 
1990), which have a radius given by 



c+2\ 



(3) 



where q is the axis ratio, and c, as defined, permits the 
ellipses to be either disky (c < 0) or boxy (c > 0). The 
shape parameters of the generalized ellipses (c, q, the cen- 
ter, and P.A.) are free parameters in the fit, but do not 
vary as a function of radius. To optimize the fit, GALFIT 
uses a Levenberg-Marquardt algorithm (Press et al. 1992), 
one that traverses down a gradient toward a minimum. 
This algorithm is very fast compared to random- walk type 
algorithms (e.g.. Simulated Annealing or Metropolis) at 
probing large parameter spaces, but has the tendency to 
"fall" into a local minimum and be content. However, rea- 
sonable initial values and probing parameter space near 
a minimum can often overcome this shortcoming. Fur- 
thermore, in light of the fact that model functions do not 
perfectly fit galaxies, the globally minimum not 
always produce the most meaningful fit. 

In order to compare with past WFPC2 and NICMOS 
studies, we have decided to use a single 2-D Nuker func- 
tion to parameterize the galaxy, even if this choice may 
not reproduce all the complex structures near the nucleus. 
Better fits can often be achieved with two or more subcom- 
ponents (L. C. Ho et al. 2001, in preparation). We create 
a Nuker model (equation 1) over the dimension of the en- 
tire NICMOS image. We then convolve it with a Tiny 
Tim PSF by following the convolution theorem. The fast 
Fourier transforms of the model image and PSF image are 
multiplied in complex space and then inverse transformed 
to obtain the final convolved model image. The generated 
model is a single bulge component described by the Nuker 
law with optimized shape parameters that provide the best 
global fit. The Nuker parameters a, (5, and 7 are con- 
strained to be positive. The effects of PSF undersampling 
in the NICMOS images can be very significant, especially 
when the profile is cuspy near the nucleus. It is important 
that the contributions from the cusp center and the sur- 
rounding portions over the pixels are extracted correctly 
(Rest et al. 2001). In GALFIT, the inner five pixels are 
divided into elliptical-polar grids centered on the galaxy, 
with radial spacings that increase in sampling closer and 
closer to the center. Thus, the integration is performed 
with increased subsampling for the inner pixels. 

The necessity of including an additional component for 
a central compact source is obvious from the surface- 
brightness profiles in most cases, consistent with the fact 
that a large fraction of the sample is spectroscopically clas- 
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sified as active galactic nuclei (AGN). In these cases wc 
added a "point" source on top of the galaxy. The point 
source is approximated by a narrow Gaussian convolved 
with the PSF. Using a Gaussian, we can determine if a cen- 
tral source is truly unresolved (intrinsic Gaussian FWHM 
< 0.5 pixel) or simply very compact (FWHM > 1 pixel). 
In cases where a point source is not evident in the surface- 
brightness profile, we determined upper limits to the point- 
source magnitudes that would have resulted in la detec- 
tion. Wc follow the use of constant boimdarics (Press 
et al. 1992) in order to set the upper limits. Adopting the 
usual convention, is defined as the sum over all pixels 
of the squares of the deviations between data and model 
image, weighted by the uncertainties. We first obtained 
the best-fit model which corresponds to the minimum 
of the fit for the galaxy without a point source. The fitting 
was then redone by introducing point sources of fixed mag- 
nitude and FWHM, but allowing all the galaxy parameters 
to vary. This results in convergence back to the best-fit 
model for the galaxy, but with increased x^ for the fit. Wc 
adjusted the point-source magnitude in each trial until the 
X^ increased by an amount equal to the reduced x^ of the 
best fit. One of the main concerns when using this method 
to obtain upper limits is that the data points are not inde- 
pendent, due to the influence of the PSF. We performed a 
simple test using PSFs of different widths and found that 
the errors on the estimated magnitudes are less than < 0.2 
magnitudes. The upper limits on the point-source magni- 
tudes, along with all other fitted parameters, are given in 
Table 2. As expected, the limits are brighter for power-law 
galaxies compared to core galaxies. 

The presence of nuclear dust lanes can lower the surface 
brightness ovvv a small range of radii close to the nucleus, 
often mimicking the presence of a point source in the sur- 
face brightness profile (eg; NGC 4150, 4261, 4374, 5273, 
and 7743). In such cases, we made best attempts to cor- 
rect for the effects of dust. The fit was first performed on 
the observed images without masking the dust features. 
The residual images were then used to locate the dust fea- 
tures, and masks were generated. The fit was subsequently 
repeated using the dust masks. If the best-fit inner slope 
cannot account for the excess light in the center, we in- 
clude a central point source. The effect of dust on the 
surface-brightness profile is very prominent in the case of 
NGC 4150, and careful masking was essential to remove 
the artificial nucleation caused by the strong, nuclear dust 
lane. 

We evaluate the quality of the 2-D fit in two ways. 
First, we perform isophote fits on the observed image and 
compare them with the 2-D fits. This is illustrated in 
the surfac;c-brightness profiles of Figure 1, where the solid 
points show the profile derived from the observed image, 
the dotted curve corresponds to the best-fitting 2-D Nuker 
profile for the bulge, and the solid curve represents the 
bulge profile plus an optional additional point source. It is 
evident that the 2-D modeling by GALFIT reproduces the 
observed azimuthally averaged profile very well, in nearly 
all cases. The second, more straightforward way to judge 
the quality of the fit is to look directly at the residual im- 
age, formed by subtracting the original image from the 2-D 
model. This is shown in the right grayscale panel of each 
galaxy in Figure 1. The residual image very effectively 



highlights the shortcoming of a single-component model. 
It can be readily seen that the Nuker function is not an 
adequate representation of the central profiles in some ob- 
jects, particularly those with strong stellar disks. For NGC 
3593 and NGC 4111, for instance, the single Nuker fit is a 
gross oversimplification given their complex morphologies. 
In NGC 3593 (Fig. If) the fitting is affected by the pres- 
ence of chaotic star-forming regions and associated dust, 
while the surface-brightness profile of NGC 4111 (Fig. Ig) 
has a number of very disky and elliptical components. In 
situations where the central point source is strong or the 
bulge profile is very peaked, Tiny Tim models the PSF 
cores well, but it may not reproduce all the fine struc- 
tures of the diffraction rings and spikes (see, e.g., NGC 
404, Fig. la; NGC 5273, Fig. Im,; NGC 5548, Fig. In). 
The residual image further accentuates faint, high-spatial 
frequency features such as dust lanes and star clusters. 

As a test of our method, we applied GALFIT to the 
F555W WFPC2 image of NGC 221 (M32) and the F547M 
image of NGC 3379, both of which have published 1-D 
Nuker fits. We performed the Nuker fit using our 2-D 
modeling procedure for the region within 1" of the galaxy 
center for M32, similar to the region used by Lauer et al. 
(1998). The fitted parameters {a = 2.13, f3 ^ 1.47, 7 = 
0.47, lib = 12.93, and rb = 0.47) are in good agreement 
with the published 1-D results (see Table 3). We also per- 
formed the fit over a larger radius (r ^10"), similar to 
the region used for our NICMOS images. The results (a 
= 4.13, P = 1.24, 7 = 0.51, lib = 12.70, and n = 0.38) 
closely agree with the parameters derived using NICMOS 
images (Table 2). In the case of NGC 3379, the Nuker fit 
for the region within the central 10" radius on the optical 
image yields a = 1.74, f3 = 1.44, 7 = 0.2L = 16.11, 
and rb = 1.86, which again agree with the published fit 
parameters. Prom these and other similar tests, we con- 
clude that the fit parameters do vary depending on the 
region chosen for the fit, but that for any given region, 
the GALFIT results closely match the parameters derived 
using 1-D fits. 

3.3. Estimation of Photometric Accuracy 

We performed aperture photometry in order to compare 
NICMOS magnitudes with ground-based NIR measure- 
ments compiled by de Vaucoulcurs & Longo (1988). In 
most cases the ground-based photometry was done with 
relatively large apertures which are beyond the area cov- 
ered by the NICMOS images, making a direct comparison 
of aperture magnitudes difficult. But for 13 galaxies, com- 
parison is possible (aperture diameter w 10"), and the 
results are shown in Figure 2. The average difference be- 
tween the two sets of data is (Amn) = {mH{HST) — 
m// (ground)) = —0.006 ±0.11 mag. This level of agree- 
ment is consistent with the assessment of Stephens et al. 
(2000), who found (Amff) = 0.048 ± 0.16 mag. 

4. MORPHOLOGY OF THE INNER REGIONS 

4.1. Dust Features 

Dust in early-type galaxies provides vital clues to un- 
derstanding their evolutionary history. Although ellipti- 
cal galaxies were once believed to be dust-free systems on 
the basis of the Hubble classification scheme, it was real- 
ized from ground-based studies that some of these systems 
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show evidence for dust (e.g., Hawardcn et al. 1981; Sadler 
& Gerhard 1985). The material revealed in the dust lanes 
may be internally generated through stellar mass loss, or 
its origin may be external, such as accretion from inter- 
actions (Forbes 1991; Knapp, Gunn, & Wynn- Williams 
1992). Nuclear dust features found within the central 
few hundred parsecs are particularly interesting from the 
standpoint of their dynamical state and their possible rela- 
tion with central activity. Because of the short dynamical 
timcscalcs in the center (^10* yr), the dust should settle 
in a plane of the galaxy in which stable orbits are allowed. 
On the other hand, if the dust is not settled, it is likely 
to have been acquired recently from an external source. 
Nuclear dust lanes are routinely detected in HST optical 
images of early-type galaxies (Laucr ct al. 1995; Forbes ct 
al. 1995; van Dokkum & Franx 1995; Verdoes Kleijn et al. 
1999; de Koff et al. 2000; Tomita et al. 2000; Tran et al. 
2001). The y-band WFPC images of early-type galaxies 
analyzed by Lauer et al. (1995) revealed dust features in 
28% of the sample, while 43% of galaxies in the study by 
Rest et al. (2001) showed dust disks and filaments. Inter- 
estingly, Rest et al. (2001) find that nuclear dust disks, 
when present, reside mostly in power-law galaxies. The 
frequency of occurrence of dust lanes appears to be con- 
nected with AGN activity. Van Dokkum & Franx (1995) 
and Verdoes Kleijn et al. (1999) find that at least 75% of 
radio-loud elliptical galaxies contain dust features in their 
F-band images. Analysis of WFPC2 images of SCR radio 
galaxies revealed dust features in one-third of the sample, 
with a correlation between the morphological distribution 
of dust and the Fanaroff & Riley (FR; 1974) classifica- 
tion (de Koff et al. 2000). When present, dust in FR I 
sources is distributed in small-scale (^2.5 kpc). sharply 
defined disks whose major axes often lie perpendicular to 
the radio jets. FR II galaxies, by contrast, display dust 
structures with a variety of sizes and morphologies. 

Early-type galaxies appear largely smooth in the NIR. 
Only ^^25% of the galaxies in our sample show evidence 
for dust in the H band. Some galaxies known to be dusty 
in i/S'T optical images (e.g., NGC 4278, 4589, and 7626) 
display smooth morphologies in the NICMOS images, il- 
lustrating that the NIR images are relatively dust free and 
can provide more accurate surface-brightness profiles. For 
instance, the profile of NGC 524 in the optical is signif- 
icantly affected by dust (Kormendy 1985a; Lauer et al. 
1995); its dust lanes form concentric shells and are clearly 
visible in the optical images of Laucr et al. (1995). By 
contrast, they are barely visible in the residual image in 
Figure lb. On the other hand, NGC 5838 has a very 
prominent dust lane that encircles the nucleus and causes 
severe obscuration of the inner regions even in the NIR. 
The WFPC2 images of this galaxy reveal another concen- 
tric dust ring located exterior to that seen in Figure In. 
Nuclear dust features are also seen in the case of NGC 
1052, 4150, and 4374, all of which host emission-line nu- 
clei at the center. NGC 3593 shows considerable amounts 
of dust and star formation activity and exhibits a chaotic 
morphology in its center. 

Since the NIR images are less sensitive to dust, explor- 
ing the relationship between dust and AGN activity in 
our sample would have to be based on the optical images. 
Dust properties for most of our sample galaxies have been 



discussed in previous works (van Dokkum & Franx 1995; 
Tomita et al. 2000; Quillen, Bower, & Stritzinger 2000). 
For those galaxies that do not have published color maps 
or residual images, we retrieved the optical images from 
the HST archive. Only two galaxies in our sample, NGC 
4026 and NGC 4417, have not been observed using HST 
in the optical. Based on examination of the archival opti- 
cal images and other published results, we find that dust 
occurs in 60% of our sample galaxies. Almost all galaxies 
with nuclear activity have dust in the central regions, the 
only exceptions being NGC 474 and NGC 5982. None of 
the nine galaxies with pure absorption-line nuclei show ev- 
idence for dust. This reinforces the results from previous 
studies on the association between dust and AGN activity. 
It also strengthens another related, emerging trend: dust 
absorption is invariably coupled with optical line emission 
in nuclear regions (Verdoes Kleijn et al. 1999; Pogge et al. 
2000). 

4.2. Nuclear Stellar Disks 

Strong nuclear stellar disks are distinctly seen in the 
residual images of NGC 2685, 3115, 4026, 4111, and 4417. 
NGC 3384 and NGC 3900 have weak disks, as seen on the 
residual images and as reflected in the positive values. 
The presence of disks in these galaxies is not surprising 
given their SO Hubble type. Two ellipticals — M32 and 
NGC 821 — show mild disklike characteristics as well. As 
was noted by Michard & Nieto (1991) from their ground- 
based images, M32 tends to have positive B4 values within 
r « 5"; they further remark that the inner disky isophotes 
become more evident at ^^1 /im. However, the WFPC2 V 
and / images of Lauer et al. (1998) do not exhibit signifi- 
cant Bi values or any presence of an inner disk of red stars. 
The residual image of NGC 821 reveals a weak disklike 
feature, and the B^^ parameter has small positive values 
over most of the semi-major axis length (see also Lauer 
1985b). Both M32 and NGC 821 possess steeply rising 
surface- brightness profiles similar to SO galaxies. The fre- 
quency of stellar disks in our sample (21%) is higher than 
that found by Lauer et al. (1995; 13%). Rest et al. (2001) 
find embedded stellar disks in a significantly larger fraction 
of their sample (51%), which is dominated by power-law 
galaxies. 

5. SURFACE-BRIGHTNESS PROFILE TYPES 

5.1. Inner Slopes 

Following the criteria proposed by Lauer et al. (1995), 
we classify a galaxy as core type when there is a significant 
flattening in the slope of the outer power-law proflle to an 
inner slope 7 which is less than 0.3. Power-law galaxies are 
defined by 7 > 0.5. Our sample has 13 core galaxies and 
17 power-law galaxies. In our sample, power-law profiles 
are mostly associated with SO galaxies, except for M32 and 
NGC 821, which are ellipticals. Interestingly, both of these 
galaxies show evidence for disky isophotes (§ 4.2). NGC 
404 and NGC 524 are the only SO galaxies with core-type 
profiles; as discussed in § 7.1, NGC 404 deviates strongly 
from the core fundamental plane, and it is very likely to 
have been misclassified. In the case of NGC 5548, its clas- 
sification as a core galaxy is somewhat uncertain due to the 
extremely bright Seyfert nucleus, whose dominant diffrac- 
tion rings contaminate the bulge profile even beyond r i=s 
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1". Out of the seven galaxies with detected stehar disks, 
all have power-law surface-brightness profiles, except for 
NGC 4111 whose classification is highly uncertain because 
the bulge cannot be modeled by a single component. In 
order to see how the isophotal shapes relate to the profile 
type, we adopt the criteria > 0.01 for disky galaxies 
and B4 < —0.01 for boxy galaxies. All nine galaxies with 
disky isophotes have power-law profiles except NGC 7457, 
which has an intermediate inner slope. Two galaxies with 
boxy isophotes have core- type surface- brightness profiles. 
Among the 15 galaxies with pure elliptical isophotes, nine 
are core-type, three have intermediate slopes and three are 
power laws. The remaining galaxies show mixed isophotal 
shapes due to the presence of dust. 

Fabcr et al. (1997) found that there is a clear dichotomy 
in the profile types with the luminosities of the galaxies. 
Core galaxies are luminous systems with My ^—20.5 mag; 
power-law galaxies are fainter and extend up to My > 
—22.0 mag. Both profile types overlap in the magnitude 
range —20.5 <,Mv ^—22.0. Within the range of resolvable 
break radii occupied by core and power-law galaxies, there 
was a clear absence of galaxies with inner-slope values from 
0.3 to 0.5. The galaxies occupied distinct regions when 7 
was plotted against log r;,, with a clear gap in the above 
range. The core galaxies in our sample are luminous sys- 
tems with boxy or pure elliptical isophotes, while power- 
law galaxies are less luminous and have disky isophotes, 
in agreement with the results of Faber et al. (1997). 

The addition of our sample to that of Faber et al. some- 
what blurs the sharp boundary between core and power- 
law galaxies (Fig. 3). Although the distribution of 7 still 
appears bimodal, four galaxies in our sample (NGC 474, 
5273, 7457, and 7626) occupy the gap region^. Rest ct 
al. (2001) provide Nuker fits of the major-axis and minor- 
axis profiles of 57 early-type galaxies and have shown that 
the cusp slopes along the two axes are similar. Including 
the inner slopes for the major axis fits from the sample 
of Rest et al. (2001) further emphasizes the continuous 
distribution of inner slopes. 

Are there any systematic errors in our profile fitting 
which could have caused some of the galaxies to fall in 
the gap region? Three of the galaxies (NGC 474, 5273, 
and 7457) have relatively strong, pointlike nuclei, and it 
is possible that the nuclei may have affected the determi- 
nation of the inner-profile slope. We do not believe this 
to be the case. Our tests indicate that the typical errors^ 
on 7 are only ±0.03. The only exception is NGC 5273, 
for which an error on 7 of ±0.20 is possible. NGC 7626 
itself has no pointlike nucleus, with an upper brightness 
limit of m^'^ ^19 mag (Table 2). Another possibility is 
that the intrinsic profiles of these objects are truly of the 
power-law variety (7 > 0.5) which, because of exceptional 
amounts of extinction, even at 1.6 /xm, have been artifi- 
cially depressed to shallower values, thereby shifting them 
into the gap region. Again, except for NGC 5273, there 
is no evidence that this is so. NGC 7626 does contain a 
kinematically distinct core (it was included in the study 
of CaroUo et al. 1997), but this kinematic attribute is ir- 



relevant to the fidelity of our profile fitting. CaroUo et al. 
(1997) obtained 7 « for NGC 7626 and classified it as a 
core type. Their flat inner slope, however, was very likely 
caused by the dust lane present in the optical image (see 
also § 5.2). Quillen et al. (2000) analyzed the H-band 
image and derived 7 = 0.46. Our analysis of the same 
data yielded a slightly shallower slope of 7 = 0.36. 

5.2. Comparison with Previous Studies 

It is of interest to compare our results with published 
material, since our 2-D fitting method differs fundamen- 
tally from methods employed in previous studies, and since 
the majority of the existing work is based on imaging at op- 
tical wavelengths. Table 3 compares the Nuker-law param- 
eters for a subset of 13 galaxies from the present work with 
those obtained from published studies based on WFPC 
and WFPC2 V^-band images. The inner-cusp slopes 7 and 
the break radii rf, are in good agreement for most galaxies 
(Fig. 4). This indicates that, to first order, there are no 
gross systematic difi'erences between our results and those 
of others. It also suggests that there are no strong color 
gradients between V and H interior to the break radius. 
By contrast, the majority of the outer slopes (3 determined 
from the NIR images tend to be steeper than the optical 
values. As the outer slopes are largely insensitive to the 
detailed treatment of the PSF, this result must reflect a 
genuine difference in the profiles between the two bands, 
one that can most plausibly be attributed to a positive 
V — H gradient toward smaller radii. 

A few cases deserve special attention. NGC 524, 4589, 
and 7626 (shown as filled circles in Fig. 4) appear consis- 
tently as outliers, most notably in the 7 and (3 plots. The 
discrepancies between the optical and NIR can be ascribed 
readily to the effects of dust, since all three galaxies show 
clear evidence for dust lanes in the optical images (Lauer 
et al. 1995; CaroUo et al. 1997). Extinction causes the 
inner slope to appear relatively flat in the optical, such 
that 7(opt) « for all three cases, but 7(NIR) system- 
atically exceeds 7(opt), by quite a signiflcant margin in 
NGC 7626. The prominent stellar disk in the SO galaxy 
NGC 3115 makes the modeling using the Nuker function 
difficult; its break radius is 85% larger in the optical than 
in the NIR. The differences in the Nuker parameters for 
M32 between our work and that of Lauer et al. (1998) is 
mainly due to the choice of the fitting radius, as discussed 
in § 3.2. 

Quillen et al. (2000) investigated the structural param- 
eters of a sample of 27 elliptical galaxies observed with 
NICMOS in the H band, as in our study. Their method of 
analysis differs from ours. They fitted a Nuker law to the 
1-D surface-brightness profiles, which were derived from 
deconvolved images. No effort was made to account for 
pointlike nuclei. Twelve of the galaxies overlap with our 
sample, and for these we find significant systematic dif- 
ferences in the derived Nuker-law parameters. The aver- 
age difference (and standard deviation) between our val- 
ues and those of Quillen et al. are as follows: (A/ih) = 
0.64 ±0.92, (Arb) = 0.071 ± 0.33, (Aa) = -0.89 ± 2.00, 

*Using early HST observations obtained with WFPC, Lauer et al. (1991) concluded that the central profile of NGC 7457 was consistent 
with a 7 Si —1.0 power law with an additional superposed pointlike nucleus. However, Lauer ct al. were only able to model the strongly 
aberrated PSF in a preliminary fashion, and we do not consider their result to be in serious conflict with ours. 

®The uncertainties represent the 68% confidence intervals, which we estimated from the use of constant boundaries (see § 3.2). 
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(A/3) = 0.14 ± 0.18, and (A7) = -0.098 ± 0.073. The 
systematic discrepancies are most simply explained by the 
treatment of nucleated galaxies. 

6. CENTRAL COMPACT SOURCES AND AGNS 

Unresolved central nuclei are present in 14 out of the 
33 sample galaxies. Excluding NGC 3593, 4111, and 5548, 
whose profile types are uncertain, we find that the fraction 
of galaxies containing unresolved nuclei is roughly equal 
between core (38%:) and power-law systems (35%). What 
is the nature of these unresolved sources? Are they AGNs 
or compact nuclear star clusters? How does their occur- 
rence relate to the large-scale properties of the host galax- 
ies? The current statistics are too meager to address these 
issues meaningfully, and we merely point out a few note- 
worthy trends. Eleven of the 14 nucleated sources arc spcc- 
troscopically classified as AGNs or closely related objects 
(eight LINERs and three Seyferts), and nearly all (5/6) 
of the "type 1" nuclei (those with visible broad emission 
lines) show evidence for point sources. Of the remaining 
objects, two have a pure absorption-line spectrum, and 
one is an H II nucleus. The majority of the nuclei, there- 
fore, appear to be associated with active galaxies, and it 
is reasonable to postulate that the 1.6 /xm emission may 
be nonstellar in origin. Lauer et al. (1995) find compact 
nuclei in 35% of their sample; they occur preferentially in 
power-law galaxies and are not associated with nonstellar 
activity. The only two cases where the central nucleus is 
an AGN turn out to be core-type galaxies. Rest et al. 
(2001) adopt more conservative criteria for identifying nu- 
clei in their sample to avoid false detections arising from 
artifacts caused by deconvolution and the effects of dust. 
They find a much lower detection rate of 13%; only two 
out of the nine nucleated galaxies have core-type profiles, 
both of which show evidence for AGN activity. 

Our statistics on the incidence of compact nuclei are 
affected in part by resolution effects and should be con- 
sidcjrecl lower limits. A handful of our galaxies were ob- 
served with the NIC3 camera, whose pixel scale is 2.5 times 
coarser than that of NIC2. The ability to discern very 
faint, photometrically distinct nuclei depends critically on 
image resolution. This is most dramatically illustrated in 
the case of NGC 4278, for which data are available from 
both cameras. The NIC3 image shows no clear indication 
of a central point source, but in the higher resolution NIC2 
image the compact nucleus emerges unambiguously from 
the underlying galaxy profile. Similarly, the NIC3 image 
of NGC 3115 used in the present study shows no evidence 
for the central star cluster found in WFPC2 images by 
Kormendy et al. (1996). 

It is of interest to note an apparent dependence of the 
incidence of nuclear ionized gas on profile type. Nearly 
all of the core-type galaxies in our sample (12/15 or 80%) 
exhibit detectable optical line emission within the central 
2"x4" (few hundred pc) above a spectroscopic equivalent- 
width limit of ~0.25 A (Ho et al. 1997a). In accordance 
with known trends of spectral class with Hubble type and 
galaxy luminosity (Ho et al. 1997b), the nuclei are classi- 
fied as AGNs or closely related objects, the majority of the 
LINER variety. Of the 17 galaxies with power-law profiles 

^"For consistency with our sample, the absolute magnitudes plotted 
couleurs et al. (1991) using distances from TuUy (1988) when available, e 
Extragalactic Database (NED), using Hq = 75 km s~^ Mpc~^. 



(including the few in the gap region), only 11 (65%) have 
detectable nebular emission, all of which are also classi- 
fied as AGNs. The marginally higher detection frequency 
of nuclear line emission among core-type systems comes 
somewhat as a surprise, given that the central regions of 
ellipticals in general emit weaker optical line emission and 
contain a smaller mass of ionized (~ 10^ K) gas than do 
lenticulars (Ho 1999b). The apparent conflict may simply 
reflect small-number statistics. 

7. DISCUSSION 

7.1. Fundamental- Plane Relations for Galaxy Cores 

Elliptical galaxies and bulges obey global correlations 
between size (rp), surface brightness (fie), velocity dis- 
persion (iJe), and luminosity - the "fundamental plane" 
(Faber et al. 1987; Dressier et al. 1987; Djorgovski & 
Davis 1987; Bender, Burstein, & Faber 1992). Luminous 
galaxies have larger effective radii (lower central densities), 
lower surface brightnesses, and higher velocity dispersions 
(Kormendy 1985b). These scaling relations have impor- 
tant implications for theories of galaxy formation, which 
must successfully reproduce them. Faber et al. (1997) 
have discussed the existence of a fundamental plane in 
the (rfc, /Ufc, c7o)-space for galaxy cores, analogous to the 
one known in the (fg, i^e, cre)-space on larger scales. In 
Figures 5 and 6, we reproduce the fundamental-plane re- 
lations for the core galaxies in our sample. The break 
radii obtained for power-law galaxies are not physically 
meaningful, and hence are not considered here, although 
these galaxies may have genuine cores on scales smaller 
than we can resolve in these observations. The log rf,- 
Mg^ (Fig. 5c) and log rt <jq (Fig. 5d) plots incorporate 
additional, nonover lapping core galaxies from the sample 
of Faber et al. (1997) and Rest et al. (2001)^°. As demon- 
strated in § 5.2, our NICMOS-based measurements of rt 
agree well with values derived from WFPC and WFPC2 
images. 

As Faber et al. (1997) have found in the optical, our 
study shows that the core parameters of elliptical and SO 
galaxies similarly follow the fundamental-plane relations in 
the NIR, again by direct analogy to well established scal- 
ing relations on global scales (e.g., Pahre, de Carvalho, 
&: Djorgovski 1998; Scodeggio et al. 1998; Mobasher et 
al. 1999). The surface brightness at the break radius and 
the break radius correlate well with the central velocity 
dispersion (Fig. 5&, 5(i), but having greater scatter with 
total galaxy luminosity (Fig. 5a, 5c) — more luminous el- 
lipticals with higher velocity dispersions have larger cores 
and lower surface brightnesses. A tight correlation exists 
between fit and ri, (Fig. 6), similar to the /ig r^ relation 
defined at the effective radius. Excluding the two outliers, 
a least-squares fit gives Hb oc ^l'^^^^'^^, with an rms scat- 
ter of 0.24 about this fit. 

Two galaxies — NGC 404 and NGC 4636 — deviate 
from the norm by their unusually low surface brightness. 
NGC 404 is a nearby {D = 2.4 Mpc), low-luminosity 
(Mg^ w —16 mag) SO galaxy which contains a very 
prominent nucleus (m™*^ = 13.5 mag; M^"^ = —13.4 
mag). The nucleus appears as a bright ultraviolet source 

in Figures 3 and 5 were derived from magnitudes given in de Vau- 
,nd otherwise from heliocentric radial velocities listed in the NASA/IPAC 
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with signatures of young, massive stars (Maoz et al. 1998). 
With 7 = 0.28 and a small break radius (r^ = 0'.'40), how- 
ever, the status of NGC 404 as a core galaxy is rather 
shaky. The inner slope depends critically on the decom- 
position of the bright nucleus, and equally acceptable fits 
can be admitted with 7 « 0.28 0.50. The classification of 
NGC 4636 as a core galaxy, on the other hand, is quite se- 
cure (7 = 0.13; Tb = 3'.' 44). The central light distribution 
is not affected by a nuclear point source (to™'^ > 23 mag) , 
and no sign of severe absorption is evident in the residual 
image (Fig. Im). The fitted Nuker parameters, moreover, 
show excellent agreement with those obtained by Faber et 
al. (1997) based on WFPC2 data (Table 3), further sug- 
gesting that dust extinction has a minimal impact on the 
photometric parameters. Thus, the displacement of NGC 
4636 from the core scaling relations appears to be genuine. 

7.2. Distribution 0/7; Bimodal or Continuous? 

Faber et al. (1997) emphasized the apparent gap in the 
distribution of the inner slopes, with no galaxies having 
0.3 < 7 < 0.5. The distribution of 7 appeared bimodal. 
Gebhardt et al. (1996) deprojected the surface-brightness 
profiles of Lauer et al. (1995) to examine the luminosity- 
density profiles, and they showed that the inner slopes of 
the latter also exhibited a bimodal distribution. However, 
the recent analysis of Rest et al. (2001) reveals that al- 
most 10% of their galaxies have asymptotic inner slopes in 
the range 0.3 < 7 < 0.5. Thus, whether the distribution 
of 7 is bimodal or not remains somewhat controversial. 

In the present study, four galaxies have inner slopes in 
the "gap region" (Fig. 3); these include the SO galaxies 
NGC 474, 5273, and 7457, and the elhptical NGC 7626. 
The presence of a bright nucleus or dust lanes can affect 
the determination of the inner slope. All three SO galax- 
ies show evidence for bright nuclei. But as mentioned in 
§5.1, our tests indicate that the only case where the point 
source may be problematic is NGC 5273, whose inner slope 
has an unusually large uncertainty (7 = 0.37±0.20). The 
IfT errors on 7 are less than 0.03 for the remaining three 
galaxies. Similarly, apart from NGC 5273, none of the 
other three galaxies show significant dust in the H band 
that could have affected the determination of 7. 

Faber et al. (1997) and Rest et al. (2001) have dis- 
cussed the ambiguity introduced in the classification of 
profile types due to distance effects. An intrinsic core 
galaxy placed at a large distance is likely to be classi- 
fied as a power-law galaxy if the break radius is close to 
the resolution limit. Similarly, for small values of a, the 
slope changes gradually with radius, and the asymptotic 
value of 7 may be reached only at radii much smaller than 
the resolution limit. Rest et al. (2001) introduced a new 
parameter, 7', which is the gradient at O'.'l derived from 
the best Nuker fit, and their classification of the surface- 
brightness profiles is based on this parameter. Faber et 
al. (1997) illustrated using M31 how distance effects can 
produce misleading 7 values. M31 is classified formally 
as a core galaxy, but it would be classified as a power- 
law galaxy if it were at the distance of the Virgo cluster. 
Intermediate values of 7, therefore, could result from mea- 
suring the inner slope too close to the break radius, es- 
pecially when a is small. In the case of the four galaxies 
in our sample with intermediate 7, the sharpness of the 



transition from the inner to outer slopes as indicated by a 
is significant, and their break radii are, for the most part, 
well resolved (rft « 0'.'33-l'.'47). Thus, it does not appear 
that the intermediate inner slopes of these galaxies can be 
attributed readily to distance effects. 

Although the introduction of the four intermediate ob- 
jects is insufficient to erase the bimodality of the 7 distri- 
bution, it is of interest to ask why previous studies have 
failed to find such objects. The studies by Lauer et al. 
(1995) and Faber et al. (1997) explicitly avoided galaxies 
known to be dusty. The selection criteria of our study, 
on the other hand, are more general, and do not bias 
against objects with dust. Since it is reasonable to sup- 
pose that cold gas content scales with dust content, one 
wonders whether our less rcstric;tive scilciction criteria may 
have uncovered certain galaxies in unusual evolutionary 
states. Although Rest et al. (2001) adopt a different def- 
inition for the inner slope, it is intriguing that they too 
find a handful of galaxies in the gap region. 

NGC 7626 contains a kinematically distinct core, a pos- 
sible relic of a former galaxy interaction or merger event. 
According to Carollo et al. (1997), however, this class 
of ellipticals does not exhibit overtly different photomet- 
ric properties on HST scales compared to kinematically 
normal ellipticals. Carollo et al. also questioned the re- 
ality of the supposed dichotomy between core and power- 
law galaxies; they suggested that the trends between in- 
ner cusp slope and global galaxy properties may instead 
be continuous. Their arguments, however, were based on 
the usage of the average logarithmic slopes of the nuclear 
profiles between 10 and 50 pc ((7phys)) instead of on the 
asymptotic inner slope of the Nuker function (7), with 
which the dichotomy between core and power-law galax- 
ies was originally proposed, and is the one adopted in this 
study. 

The remaining three sources with intermediate 7 are 
SO galaxies. NGC 5273 is uncertain, as described above. 
However, apart from being nucleated, itself not an uncom- 
mon attribute, neither NGC 474 nor NGC 7457 show any 
particularly noteworthy characteristics. They are not ex- 
ceptionally dusty; their global optical colors appear typical 
of SO galaxies; and their central stellar velocity dispersions 
roughly conform to the Faber- Jackson (1976) relation for 
SOs, indicating normal mass-to-light ratios. 

The present statistics on galaxies with intermediate cusp 
slopes are clearly too meager to warrant excessive specula- 
tion on their physical nature. Larger samples are needed 
to place these objects in the context of current forma- 
tion scenarios for power-law and core galaxies (Faber et 
al. 1997). 

7.3. Isothermal Cores 

A remarkable finding that emerged from high-resolution 
ground-based studies is the virtual absence of isothermal 
cores in "normal" elliptical galaxies (Lauer 1985a; Kor- 
mendy 1985a, 1985b). The handful of objects with can- 
didate isothermal cores all turn out to be bright cluster 
galaxies. HST studies confirmed that King models with 
constant-density cores provide a poor description of the 
observed central surface-brightness profiles of ellipticals. 
Galaxies which showed distinct cores had inner profiles 
with a central rising cusp. Lauer et al. (1995) empha- 
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sised that even galaxies with very shallow cusps in the 
inner region of surface-brightness profile have steeply ris- 
ing central luminosity-density profiles, thereby ruling out 
the possibility of isothermal cores. 

Six galaxies in our sample formally have 7 < 0.05; the 
errors on 7 are small, typically less than ±0.02. Both NGC 
4261 and NGC 4278 have unresolved nuclei, making their 
inner slopes less reUable. NGC 524 and NGC 4472 have 
weak central cusps, with 7 = 0.03 and 0.04, respectively. 
(Ferrarese et al. 1994 fitted the optical HST profile of 
NGC 4472 with an isothermal model, but the model does 
not fit the observed profile well, and dust is clearly present 
in their image.) Excluding these cases, the two remaining 
objects, NGC 4291 and NGC 4406, are excellent candi- 
dates for galaxies with isothermal cores. The core of NGC 
4406 is very well resolved {r^ = I'.'OO), and the inner slope 
is identically zero (Icr error ±0.003). The image appears 
very smooth, showing no evidence of a distinct nucleus, 
strong dust absorption, or other irregularities. Kormendy 
(1985a) already demonstrated that the central region of 
NGC 4406 is well fitted by an isothermal profile; beyond 
~6" the isothermal model is inadequate. Our analysis, 
based on observations at a redder bandpass and much 
higher resolution, lends greater confidence to Kormendy's 
finding. The core in NGC 4291 is slightly more compact 
(fb = 0'.'48), but there is also no reason to suspect that the 
surface-brightness profile might be anomalous, and it has 7 
= 0.02±0.02. The ease for an isothermal core in NGC 4291 
is thus somewhat less convincing than in NGC 4406, but 
we consider it a possible candidate nonetheless. Within 
the errors, the logarithmic slope of the luminosity-density 
profile is nearly zero in both these galaxies. Neither object 
is a bright cluster galaxy. 

A few galaxies in the Faber et al. (1997) sample have 
flat, central surface-brightness profiles with < 7 < 0.02 
(Fig. 3). Only two of them have cores with constant 
luminosity-density, namely NGC 1600 and NGC 4889 
(Gebhardt et al. 1996). NGC 4889 is one of the brightest 
cluster galaxies in Coma, but NGC 1600 is not a bright 
cluster galaxy. Thus, while isothermal cores are generally 
thought to occur exclusively in bright cluster galaxies, a 
minority of objects appear to violate this trend. 

7.4. Properties of Central Unresolved Sources 

Until HST images became available it was difficult to 
identify and quantify photometrically distinct nuclei em- 
bedded in the bulges of nearby galaxies. The difficulty 
of obtaining bulge-free magnitudes for AGNs is clearly 
borne out by the lack of reliable luminosity functions for 
nearby AGNs (Krolik 1998). Two extreme examples in our 
sample illustrate the point. The luminous elliptical NGC 
4278 contains a weak {m^'^ = 17.2 mag), intrinsically 
faint (Af™*^ = —12.8 mag) nucleus which is overwhelm- 
ingly swamped by host-galaxy light. By contrast, the well 
known Seyfert nucleus in NGC 5548 (m™^ = 11.9 mag; 
M^"'= = —22.2 mag) accounts for a significant fraction of 
the total luminosity. We have applied our 2-D decomposi- 
tion method self consistently to measure or constrain the 
strength of the nuclear component. Excluding NGC 5548, 
detected nuclei have magnitudes in the range to™'^ = 12.8 
to 17.4 mag, or M^" = -12.8 to -18.4 mag. Undetected 
sources have limits as low as m^'^ » 23 mag and M^"^ » 



—6.5 mag. 

The physical nature of the compact nuclei is not well 
constrained. All of the sources appear unresolved at the 
resolution of the NICMOS images, which corresponds to 
FWHM diameters of 14 and 18 pc for NIC2 and NIC3, 
respectively, for the sample median distance of 17 Mpc. 
These scales are insufficient to distinguish between a truly 
pointlike nucleus and a compact star cluster. With a few 
exceptions, we also have no information on the spectral 
properties of the sources on nuclear scales. However, as 
discussed in § 6, the association of the nuclei with AGN 
classifications (as determined from ground-based spectra) 
suggests that they may indeed be nonstellar in origin. If 
this is the case, the NIR continuum emission might scale 
with the optical line luminosity (Fig. 7a). The two vari- 
ables are marginally correlated at the 94% confidence level 
according to the generalized Kendall's r test (Isobe, Feigel- 
son, & Nelson 1986). However, the significance disappears 
altogether when we properly account for the mutual de- 
pendence of the two variables on distance; the partial 
Kendall's r test (Akritas & Siebert 1996) cannot reject 
with 52% probability the null hypothesis that there is no 
correlation. The lack of a clear relation between the NIR 
continuum and optical line emission, on the other hand, 
is hardly surprising, since the l.Q-fxm window is far re- 
moved from the ionizing part of the spectrum. Moreover, 
the spectral energy distributions of low-luminosity AGNs 
exhibit a diversity of forms (Ho 1999c), rendering any cor- 
relation between the NIR and far-ultraviolet bands non- 
trivial. 

Finally, we note that the NIR luminosity of the nuclear 
sources does not depend on the central black-hole mass. 
Recent observations suggest that massive black holes may 
be ubiquitous in bulge-dominated galaxies (Magorrian et 
al. 1998; Richstone et al. 1998; Ho 1999a) and that their 
masses arc tightly correlated with the stellar velocity dis- 
persion of the bulge (Ferrarese & Merritt 2000; Gebhardt 
et al. 2000). A plot of M™'^ versus central velocity dis- 
persion reveals no correlation (Fig. 76). To the extent 
that the NIR continuum traces emission from the accre- 
tion flow, this is indicative of the enormous range spanned 
by the radiative output of massive black holes in nearby 
early- type galaxies. 

8. SUMMARY 

We analyzed iJ-band NICMOS images of the central re- 
gions of 33 nearby early-type (E, SO, and SO/a) galaxies 
using a new 2-D fltting technique. Our method models the 
galaxy bulge or spheroid using the Nuker function and si- 
multaneously accounts for the PSF and a possible nuclear 
point source. Our main results are as follows: 

(1) Early- type galaxies have relatively smooth light dis- 
tributions in the H band, showing dust features only in 
very few cases, mostly associated with galaxies hosting 
emission-line nuclei. The optical y-band WFPC2 images 
show evidence for dust absorption in all the galaxies in the 
sample that have nuclear activity, while dust is absent in 
galaxies with no nuclear line emission. 

(2) The Nuker parameters derived from the NICMOS 
data generally show good agreement with published val- 
ues based on optical WFPC2 images. 

(3) Consistent with earlier studies, galaxies with distinct 
cores tend to be luminous systems with boxy isophotes, 
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whereas less luminous galaxies with disky isophotes tend 
to have power-law surface-brightness profiles with no 
clearly defined break. However, we do not find a clear 
dichotomy in the distribution of the inner cusp slopes of 
the Nuker function. Specifically, there arc several objects 
with 7 between 0.3 and 0.5, a region previously thought 
to be empty. 

(4) Core galaxies obey the core fundamental-plane re- 
lations in (rf,, ^i,, (To)-space, with more luminous, higher 
velocity dispersion objects having larger cores and lower 
surface brightnesses. 

(5) Two galaxies (NGC 4291 and NGC 4406) have con- 
stant central surface-brightness profiles. They are good 
candidates for having constant-density, isothermal cores. 
Neither galaxy is in a rich cluster, unlike other systems 
previously found to host candidate isothermal cores. 

(6) Approximately half of the galaxies in our sample are 
nucleated, roughly evenly split between core and power- 
law galaxies. Photometrically distinct nuclei are especially 
prevalent in "type 1" AGNs (those with detected broad 



emission lines). 

(7) The pointlike nuclei have magnitudes in the range 
m'^'' = 12.8 to 17.4 mag, which correspond to M^" = 
— 12.8 to —18.4 mag. We find no significant correlation 
between the NIR luminosities and the optical line lumi- 
nosities. The strength of the nuclear NIR emission also 
appears to be unrelated with the central stellar velocity 
dispersion, an indirect indicator of the black-hole mass. 
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APPENDIX 

NOTES ON INDIVIDUAL OBJECTS 

This Appendix gives a short description of each object. Most of our sample galaxies have been studied extensively in 
various contexts. We do not attempt to provide a detailed description of each galaxy, but instead focus on the significant 
morphological features, peculiarities, or any noteworthy difficulties encountered in the 2-D modeling. 



NGC 221 (M32). — The significantly positive values within r ~ 5" suggest that a disky component is present in 
M32, although the residual image does not show clear evidence for it. As in the optical (Lauer et al. 1998), the galaxy is 
remarkably smooth and featureless in the NIR. 

NGC 404- — Apart from the strong nuclear point source, the NICMOS image of NGC 404 is relatively featureless. The 
excess positive emission in the center of the residual image is due to a slight mismatch between the point source and the 
PSF. WFPC2 optical images (Pogge et al. 2000) show a spiral pattern extending all the way into the nucleus, vague hints 
of which can be seen in the NIR residual image. Interestingly, the best fit for the surface-brightness profile gives 7 = 0.28, 
officially in the domain of core galaxies. This is unexpected given the low luminosity and low velocity dispersion of the 
galaxy. As seen from Figures 5 and 6, NGC 404 deviates strongly from the locus of core galaxies in the core-parameter 
relations. It is likely that our estimate of 7 has been affected by the bright nucleus. 

NGC 4^4- — Although NGC 474 is known to be a shell galaxy with obvious signatures of a merger remnant at large radii 
(Turnbull, Bridges, & Carter 1999), the inner regions appear undisturbed, showing no evidence for significant structure 
or dust, in either the V ox H bands. The low-level features in the residual image are due to a slightly mismatched PSF 
for the nucleus and deviations from the Nuker-law fit. 

NGC 524- The presence of a concentric dust pattern in NGC 524 causes the inner profile to flatten considerably in 
the optical. The dust rings are less prominent in the NIR, but they can still be seen in the residual map. Our 2-D Nuker 
fit of the i7-band image gives 7=0.03±0.01. Using the same data, Quillen et al. (2000) obtained a much steeper slope of 
7 = 0.25. We performed the 2-D fit using their derived Nuker parameters and find that even though the observed profile 
is well reproduced for the range 0'.'3-5'.'0, the fit deviates considerably from the observed profile in the inner and outer 
regions. Within r = O'.'l, 7 = 0.25 makes the inner regions much brighter than observed. 

NGC 821. — No evidence for dust is seen in the archival U-band image. The features in the residual image are due 
to slight deviations from the Nuker-law fit. The high ellipticity and positive B4 values over most of the semi-major axis 
length suggest that a disk is present. 

NGC 1052. — The faint structure in the residual image is most likely associated with the gas and dust complex seen in 
optical images. Pogge et al. (2000) argue that the morphology of the optical emission is reminiscent of ionization cones 
found in Seyfert galaxies. 

NGC 2685. — This SO galaxy shows a strong disk within r « 2". As evident in the residual image, the disk is poorly 
modeled with a single-component Nuker function. The dust patches in the archival V image are largely absent in the 
NIR, except for the extended filament located ~5" to the northeast of the nucleus. 

NGC 3115. — The residual image shows a nuclear disk in addition to the large-scale disk; neither is well fitted with 
the Nuker law. NGC 3115 is known to have a F « 17 mag central star cluster in the optical (Kormendy et al. 1996). 
Our NIC3 image does not detect the cluster, presumably because of its low contrast and the poor pixel resolution. The 
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optical images show no dust features (Lauer et al. 1995; Tomita et al. 2000). 

NGC 3379. — The weak dust ring seen in optical images (van Dokkum & Franx 1995; Tomita et al. 2000) is not 
present in the NIR. The features in the residual image are due to slight deviations from the Nuker-law fit. 

NGC 3384- — The optical images show no dust features (van Dokkum & Franx 1995; Tomita et al. 2000). The residual 
image highlights the weak nuclear and large-scale disks, which are poorly fitted by the Nuker law. The pointlike nucleus 
is also incompletely removed by the PSF. 

NGC 3593. — The central region has a very chaotic morphology from the copious dust lanes and star-forming regions. 
The 2-D Nuker fit is obviously corrupted by the dust features, and we consider it to be unreliable. Fortunately, the strong 
nuclear source is seen quite distinctly, and we can measure its magnitude with little complication. 

NGC 3900. The archival V image shows patchy dust, but this is not present in the NIR image. As seen from the 
residual image and the radial variation of the S4 parameter, there is a nuclear disk within the inner O'.'S. The profile of 
NGC 3900 would probably be better modeled by multiple components. To achieve an acceptable fit with the Nuker law, 
we restricted the fit to r < 3". 

NGC 4026. — The strong disk in the center complicates the modeling. There are no obvious dust features in the 
residual image. No optical imago is available. 

NGC 4111- — The central region appears to be composed of a point source, a peanut-shaped bulge, and a dominant 
edge-on disk. Since parameterizing this multi-component light distribution by a single Nuker law is an extreme oversim- 
plification, we use the 2-D fit only as a guide to extract the brightness of the nucleus. We consider the specific values of the 
Nuker parameters to be unreliable. The archival V^-band image shows fan-shaped dust features emanating perpendicularly 
from the disk. This is barely visible in the residual NIR image. 

NGC 4143. — The archival V-band image shows patchy dust obscuration close to the nucleus, qualitatively similar to 
the low-level features in our residual image. 

NGC 4150. - A strong nuclear dust lane causes considerable obscuration of the central region, both in the optical and 
in the NIR. The dust lane was carefully masked while doing the 2-D fit. 

NGC 4261. — This galaxy is remarkable for its high degree of boxincss. The nuclear dust disk of NGC 4261 appears 
very prominently in optical iJST images (Jaffe et al. 1993; Verdoes Kleijn et al. 1999), but it is much less conspicuous in 
the NIR. A pointlike nucleus is definitely required to fit the inner light profile. Quillen et al. (2000) obtained a steeper 
inner cusp slope because they did not account for the excess emission from the nucleus. 

NGC 4278. — The irregular dust patches of the optical images (van Dokkum & Franx 1995; CaroUo et al. 1997; Tomita 
et al. 2000) are barely visible in the NIR residual map. The central point source is slightly mismatched by the PSF. 

NGC 4291. — The residual image does not show any special feature except for a weak quadrupole pattern, which 
results from the boxy outer isophotes (Ebneter, Davis, & Djorgovski 1988). No dust is obviously present in an archival 
F-band image. The surface-brightness distribution steepens with respect to the best-fit Nuker profile for r ^3"; we thus 
restricted the fit to r < 3". 

NGC 4374. — The multiple dust lanes of NGC 4374 (M84) are clearly visible in the residual image and have been 
extensively discussed in the literature (van Dokkum & Franx 1995; Bower ct al. 1997; Verdoes Kleijn et al. 1999). 
Extensive Ha emission is associated with the dust. The embedded nucleus, detected in the optical (Bower ct al. 1997) 
and in the NIR (this study), is significantly reddened by the dust (Bower et al. 1997). 

NGC 44O6. — This galaxy is well known for its photometric and kinematic peculiarities. It exhibits minor-axis 
rotation and has a kinematically distinct core (Forbes et al. 1995). Surface-brightness profiles obtained from ground- 
based observations under good seeing showed that it has an isothermal core (Kormendy 1985a). The HST profile is 
remarkably fiat, both in the optical in the NIR; our fit formally gives 7 = 0.00±0.003. Little or no color gradient has 
been measured by Carollo et al. (1997) and Tomita et al. (2000), suggesting that the flattening of the central light proflle 
is probably intrinsic and not due to dust extinction. We consider this galaxy to be a strong candidate for hosting an 
isothermal core (§ 7.3). 

NGC 441'^- — The 2-D fit is affected by the strong nuclear disk. No optical //ST images are available. 

NGC 4472. — As in NGC 4406, the surface-brightness profile of NGC 4472 is very fiat (7 = 0.04) in the center and 
the isophotal parameters show large variations. Irregular, patchy obscuration affects the center of the optical images 
(van Dokkum & Franx 1995; Tomita et al. 2000), but our NIR residual image appears quite smooth. 

NGC 4589. — The complex gas and stellar kinematics of NGC 4589 suggest that it is a merger remnant (MoUenhoff 
& Bender 1989). The dust filaments which traverse the center of the optical image (Tomita ct al. 2000; Quillen et al. 
2000) are also visible in our NIR residual map. The light distribution beyond ~5" deviates from the Nuker profile, and 
we confined our fit to the region interior to this radius. 

NGC 4G36. Van Dokkum & Franx (1995) find irregular dust lanes in the optical image, but our NICMOS image 
is extremely smooth. Although the ellipticity and ,64 parameter are both well behaved, the position angle shows large 
variations. As discussed in § 7.1, NGC 4636 deviates markedly from the fundamental-plane relations for galaxy cores; 
its surface brightness is too low compared to what is expected for galaxies of similar absolute luminosity and velocity 
dispersion. Likewise, its break radius is systematically larger than expected, implying an unusually diffuse core. 

NGC 5273. — The galaxy hosts a bright Seyfert nucleus. The dust absorption pattern seen in the archival optical 
image corresponds closely to the structures in the NIR residual map. 

NGC 5548. — The extremely bright Seyfert nucleus dominates over the bulge fight of NGC 5548 for radii ^l"-2". The 
Nuker parameters are thus quite uncertain. No useful (unsaturated) optical image was found in the archive. 
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NGC 5838. — The thick nuclear dust ring partly occults the nucleus and compromises the Nuker fit for the bulge. We 
were unable to improve the fit by masking out the dust ring. No estimate of the nuclear point source is given. A second 
concentric dust ring at larger radii is visible on archival optical images, and it can be faintly seen in the iJ-band residual 
image. 

NGC 5982. -— There is no evidence for dust or a point source in this galaxy, which is known to contain a kincmatically 
distinct core (Forbes et al. 1995). The isophotal structure, however, is very complicated. The central isophotes are 
perfectly circular, but beyond ~1" they become increasingly boxy. This leads to the quadrupole pattern in the residual 
image. 

NGC 6340. — The archival F-band image shows a faint dust lane cutting across the nucleus, but there is no evidence 
of it in the NIR. 

NGC 7457. — The nuclear compact source is prominent both in the NIR and in the optical (Lauer et al. 1991). There 
is no evidence for dust in the optical images (Tomita et al. 2000). The faint features in the NIR residual image come 
from a slight mismatch between the model and the data in the region r « O'.'5-l". 

NGC 7626. — The warped nuclear dust lane seen in the optical images is likely to be the culprit for flattening the 
optical profile toward the center (Carollo et al. 1997). The dust lane is absent from our NIR residual image, and we 
obtain a much steeper central cusp slope (7 = 0.36 instead of 7 = 0.0). 

NGC 7743. — The bright Seyfert 2 nucleus of NGC 7743 is surrounded by clumpy dust, clearly seen both in our 
residual image and in an archival l^-band image. The galaxy ccntroid is badly off-centered on the NIC2 image and lies in 
the upper right quadrant. Thus, only the region within r = 6" was used for the analysis presented here. 
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FIGURE CAPTIONS 

Fig. la, b. — (a) NGC 221 and NGC 404. (6) NGC 474 and NGC 524. Top: The observed NICMOS F160W image 
{left) and residual image {right). Positive values are dark, and negative values are white. The images are 10" x 10" for 
NIC2 data and 25" x 25" for NIC3 data, centered on the galaxy, with North oriented up and East to the left. Bottom,: 
Results of the isophotal analysis. The panels plot the variation of surface brightness /i//, cUipticity e, P.A., and the shape 
parameter B^^ along the semi-major axis. The best-fitting Nuker function, derived from the 2-D model, is overplotted for 
^iH with {solid line) and without {dotted line) a central point source. In the bottom three panels, the vertical dashed line 
indicates the region interior to which the points are strongly affected by the PSF. 

Fig. Ic, d. — (c) NGC 821 and NGC 1052. {d) NGC 2685 and NGC 3115. As in Fig. la, lb. Top: The observed 
NICMOS F160W image {left) and residual image {right). Bottom: Results of the isophotal analysis. 

Fig. le, /. — (e) NGC 3379 and NGC 3384. {f) NGC 3593 and NGC 3900. As in Fig. la, 1&. Top: The observed 
NICMOS F160W image {left) and residual image {right). Bottom: Results of the isophotal analysis. 

Fig. Ig, h. — {g) NGC 4026 and NGC 4111. {h) NGC 4143 and NGC 4150. As in Fig. la, lb. Top: The observed 
NICMOS F160W image {left) and residual image {right). Bottom: Results of the isophotal analysis. 

Fig. li, j. — {i) NGC 4261 and NGC 4278. (j) NGC 4291 and NGC 4374. As in Fig. la, lb. Top: The observed 
NICMOS F160W image {left) and residual image {right). Bottom: Results of the isophotal analysis. 

Fig. Ik, I. — {k) NGC 4406 and NGC 4417. {I) NGC 4472 and NGC 4589. As in Fig. la, 16. Top: The observed 
NICMOS F160W image {left) and residual image {right). Bottom: Results of the isophotal analysis. 

Fig. Im, n. (m) NGC 4636 and NGC 5273. (n) NGC 5548 and NGC 5838. As in Fig. la, 1&. Top: The observed 
NICMOS F160W image {left) and residual image {right). Bottom,: Results of the isophotal analysis. 

Fig. lo, p. — (o) NGC 5982 and NGC 6340. {p) NGC 7457 and NGC 7626. As in Fig. la, 16. Top: The observed 
NICMOS F160W image {left) and residual image {right). Bottom: Results of the isophotal analysis. 

Fig. Iq. — {q) NGC 7743. As in Fig. la, 16. Top: The observed NICMOS F160W image {left) and residual image 
{right). Bottom: Results of the isophotal analysis. 

Fig. 2. — Comparison of aperture photometry derived from NICMOS and ground-based images. The measurements 
have been made with ~10" diameter apertures. The solid line denotes equality. The average difference between the two 
sets of data is (Ami/) = {mH{HST) — mi/(gi'ound)) = — 0.006±0.11 mag. 

Fig. 3. — Dependence of the inner cusp slope 7 on (a) absolute magnitude of the galaxy {M^^) and (6) break radius 
(rb). The distribution of 7 is shown in panel (c). Filled symbols represent core galaxies, and open symbols denote 
power-law galaxies. The galaxies from this study are plotted as circles, those from Faber et al. (1997) are shown as 
triangles, and those from Rest et al. (2001) as squares. Objects in our sample which fall in the "gap region" (0.3 < 7 < 
0.5) appear as asterisks, while those from Rest et al. (2001) are shown by plus signs. Two galaxies which fall in this 
region from Faber et al. (1997) are shown by crosses. 

Fig. 4- — Comparison of the values of (a) 7, (6) /?, and (c) derived from NICMOS images to those derived from 
optical images. The solid line denotes equality. Filled circles correspond to NGC 524, NGC 4589 and NGC 7626, which 
are known to have dust in the central regions as seen on optical images. The central profiles of these objects most likely 
were depressed in the optical by dust extinction, leading to 7(opt) « 0. 

Fig. 5. — Central-parameter relations for core galaxies. The core parameters (rj, and Hb) are well correlated with 
the central stellar velocity dispersion (o"o) and, with greater scatter, with the total _B-band absolute magnitude {Mg^). 
Panels (a) and (6) plot only data from this study {circles), whereas panels (c) and {d) include measurements from Faber 
et al. (1997; triangles) and Rest et al. (2001; squares). The outliers NGC 404 and NGC 4636 are labeled. 

Fig. 6. — Fundamental-plane relation between jib and r(, for core galaxies. The outliers NGC 404 and NGC 4636 are 
labeled. 

Fig. 7. — Correlation of nuclear point-source magnitude with (a) extinction-corrected Ha luminosity and (6) central 
stellar velocity dispersion. Filled symbols denote nuclei in core galaxies, open symbols represent nuclei in power-law 
galaxies, and asterisks mark objects with intermediate inner slopes (0.3 < 7 < 0.5). Upper limits are indicated with 
arrows. Two objects with uncertain nuclear profiles appear as triangles. 
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TABLE 1: Properties of the Sample 



Galaxy Hubble D Rt log Ihc Tn Spectral Exposure EST GO 

Name Type (Mpr) (mag) (mag) (ergs~^) (kms~^) Class Time (s) Camera Program 

(1) (2) (.1) (4) (5) (6) (7) (8) (9) (10) (11) 



NGC 221 


E2 


0.7 


9.03 


-15 


.51 


<36.19 


77 


A 


16 


NIC2 


7171 


NGC 404 


SO-: 


2.4 


11.21 


-15. 


.98 


37.76 


54 


L2 


320 


NIC2 


7330 


NGC 474 


SO 


32.5 


12.37 


-20 


.59 


38.82 


169 


L2:: 


128 


NIC2 


7268 


NGC 524 


S0+ 


32.1 


11.30 


-21 


.36 


38.58 


242 


T2: 


160 


NIC2 


7886 


NGC 821 


E6? 


23.2 


11.67 


-20 


.11 


<38.13 


207 


A 


160 


NIC2 


7886 


NGC 1052 


E4 


17.8 


11.41 


-19 


.90 


39.45 


222 


L1.9 


160 


NIC2 


7886 


NGC 2685 


SBO+ pec 


16.2 


12.12 


-19 


.23 


.19.21 


114 


S2: 


192 


NIC 3 


7919 


NGC 3115 


SO— spin 


6.7 


9.87 


-19 


.39 


<37.80 


264 


A 


192 


NIC 3 


7919 


NGC .1379 


El 


8.1 


10.24 


-19 


..16 


37.94 


209 


L2:: 


128 


NTC2 


7453 


NGC 3384 


SBO-: 


8.1 


10.85 


-18 


.79 


<37.56 


170 


A 


128 


NIC2 


7453 


NGC 3593 


so/a: 


5.5 


11.86 


-17 


.25 


38.90 


77 


H 


192 


NIC3 


7919 


NGC 3900 


S0+ 


29.4 


12.20 


-20 


.17 


.18.09 


118 


7,2: 


128 


NIC2 


7331 


NGC 4026 


SO spin 


17.0 


11.67 


-19 


.56 


<38.10 


195 


A 


192 


NIC3 


7919 


NGC 4111 


S0+: spin 


17.0 


11.63 


-19 


.55 


.19.84 


139 


L2 


192 


NIC 3 


7919 


NGC 4143 


SABO 


17.0 


11.65 


-19 


.25 


38.69 


270 


L1.9 


320 


NIC2 


7330 


NGC 4150 


SO? 


9.7 


12.44 


-17 


.53 


38.18 


85 


T2 


160 


NIC2 


7886 


NGC 4261 


E2+ 


35.1 


11.41 


-21 


.37 


.39.82 


326 


L2 


192 


NTC2 


7868 


NGC 4278 


E1+ 


9.7 


11.09 


-18 


.96 


39.17 


250 


LI. 9 


160 


NIC2 


7886 


NGC 4291 


E 


29.4 


12.43 


-20 


.09 


<38.38 


278 


A 


160 


NIC2 


7886 


NGC 4374 


El 


16.8 


10.09 


-21 


.12 


39.31 


296 


L2 


192 


NIC2 


7868 


NGC 4406 


E3 


16.8 


9.83 


-21 


.39 


<37.82 


250 


A 


128 


NIC2 


7453 


NGC 4417 


SBO: spin 


16.8 


12.00 


-19 


.12 


<37.88 


84 


A 


192 


NIC3 


7919 


NGC 4472 


E2 


16.8 


9.37 


-21 


.80 


37.59 


303 


S2:: 


128 


NIC2 


7453 


NGC 4589 


E2 


30.0 


11.69 


-20 


.71 


39.50 


228 


L2 


160 


NIC2 


7886 


NGC 4636 


E0+ 


17.0 


10.43 


-20 


.72 


38.27 


207 


LI. 9 


160 


NIC2 


7886 


NGC 5273 


SO 


21.3 


12.44 


-19 


.26 


38.70 


52 


S1.5 


320 


NIC2 


7330 


NGC 5548 


so/a 


68.6 


13.30 


-21 


.32 


40.22 




S1.5 


64 


NIC2 


7172 


NGC 5838 


SO- 


28.5 


11.92 


-20 


.56 




290 


T2:: 


128 


NIC2 


7450 


NGC 5982 


E3 


38.7 


12.04 


-20 


.89 


38.46 


256 


L2:: 


160 


NIC2 


7886 


NGC 6340 


so/a 


22.0 


11.87 


-20 


.04 


.18.50 


137 


L2 


128 


NTC2 


7331 


NGC 7457 


SO-? 


12.3 


12.09 


-18 


.69 


<37.24 


77 


A 


128 


NIC2 


7450 


NGC 7626 


E: pec 


45.3 


12.16 


-21 


.23 


38.81 


273 


L2:: 


160 


NIC2 


7886 


NGC 7743 


SBO+ 


24.4 


12.38 


-19 


.78 


40.24 


85 


S2 


320 


NIC2 


7330 



Note. — Col. (1) Galaxy name. Col. (2) Hubble type. Col. (3) Adopted distance as given in Tully 1988 for D < 40 
Mpc, and otherwise derived from the heliocentric radial velocity and Ha = 75 km s^"'^ Mpc^"'. Col. (4) Total apparent 
B magnitude. Col. (5) Total absolute B magnitude, corrected for Galactic extinction and to face-on inclination. Col. 
(6) Extinction-corrected Inminosity for the narrow component of the Hq emission line. Col. (7) Central stellar velocity 
dispersion. Col. (8) Spectral class of the nucleus, where A — absorption-line nucleus, H — H II nucleus, L — LINER, 
S = Seyfert, T = "transition object" (LINER/H ll), 1 = type 1,2 = type 2, and a fractional number between 1 and 
2 denotes various intermediate types; uncertain and highly uncertain classifications are followed by a single and double 
colon, respectively. Col. (9) Exposure time. Col. (10) NICMOS camera. Col. (11) GO program number. Data for 
cols. (2)-(6) and (8) taken from Ho et al. 1997a. Some of the Ho luminosities in col. (6) have been updated with more 
accurate values than those given in Ho et al. 1997a. Data for col. (7) taken from McElroy 1995, except for NGC 4143, 
which comes from Di Nella et al. 1995. Nelson & Whittle 1995 give a velocity dispersion for NGC 5548, but it is highly 
uncertain. 
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TABLE 2: Fitted Parameters 



Galaxy 


Profile 


y"6 


n 


a 


/3 


7 




Notes 


Name 


Type 


(mag arcsec"^) 


(arcsec) 








(mag) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 221 




9 49 


\j . •)•) 


4 CP, 


1 26 


n 

\) .f.j\> 


^1 7 n 

^ -L i -U 






n 
1 


1 9 1 f! 

-LO . -LU 


n 40 


u . uo 


1 ^18 

-L . 00 


n 9^ 

u .zo 


1 9 ^19 

-LO . OZ 




NHP 474 


\ 
/ 


14 '^9 
1^ . oz 


1 47 


1 9'^ 
1 . zo 


1 Qfl 
1 . y u 


U -O i 


lu . uo 


a 




1 1 


lO . i o 


1 41 
1.41 


u . uo 


1 RQ 
I . ( ) y 


U - UO 


^ J- y . 




NHP 891 




1 9n 

lO . zu 


1 1 9 

1 . 1 z 


1 nn 

1 . uu 


I . oy 


n R4 

U .U4 


^ Z VJ .u 




NCiC in'i9 


n 
1 1 


11.91 


n 40 


1 n^i 

-L . uo 


1 .43 


n 1 1 


1 5.49 




NCiC 9fiSfS 




14 18 


z . oo 


1 f!Q 
-L . u y 


1 52 


n 7*^ 

u . 1 


^1 8 n 

^ -L -U 




NHP SI 1 




1 S 

-LO . oo 


u . o y 


1 1 S 

1 . lO 


1 sn 

-L . OU 


n 7*^ 

u . 1 


^1 

^ -L .0 




i\ vj\^ oo / y 


1 


1 9 sn 

LZ . OU 


1 • uo 


1 89 
I . o z 


1 4^1 


n 1^ 

u • 1 


/ z . u 




i\vj\^ ooo^ 




n 1 8 

I- O . I o 


9 4*^ 

Z .40 


o ■ 'ty) 


1 ^18 
I .00 


n R4 

U .U4 


14 79 




i\VjrV.j ooyo 














1 T 99 

L .J . ZZ 


u 


i\vj\^ oyuu 




1 9 '^'i 
iz . oo 


n 9*^ 

U • zo 


n 9Q 
u . zy 


1 fifi 
1 . UU 


n ^1 

U -U 1 


^ zu .0 




i\VJt\^ ^uzu 




12.40 


n ^4 

U -04 


n 88 
u . oo 


-L . OU 


U -UO 


^9n n 

^ zu .u 




Nnr 41 1 1 














12 79 








14 26 


1 1 

. 1 1 


1 26 


2 18 


n 

u -uy 


1 P, np 

-LU . UU 




Nnr* 41 




1 9 Qfi 
iz . yo 


u . uo 


1 9'^ 
1 . zo 


1 fi7 

1 . U i 


u .uo 


^1 Q 8 
^ 1 y .0 




NHP 49(^1 

iNvJtVJ 4:ZUl 


1 1 


lo . oo 


1 P,9 
1 . uz 


9 '^8 
z . oo 


1 4'^ 
1 .^0 


n nn 

U -UU 


1 7 '^fi 

1 i . OU 




NGC 4278 


n 


12.80 


0.97 


1.63 


1.39 


0.02 


17.15 




NGC 4291 


n 


12.3! 


0.48 


2.07 


1.48 


0.02 


>23.0 




NGC 4374 


n 


13.35 


2.39 


2.15 


1.50 


0.13 


16.17 




NGC 4406 


n 


12.98 


1.00 


3.31 


1.16 


0.00 


>23.0 




NGC 4417 


\ 


14.74 


3.66 


0.87 


1.77 


0.71 


>17.0 




NGC 4472 


n 


13.52 


2.63 


1.89 


1.29 


0.04 


>23.0 




NGC 4589 


n 


12.05 


0.21 


1.09 


1.18 


0.11 


>23.0 




NGC 4636 


n 


14.60 


3.44 


1.69 


1.56 


0.13 


>23.0 




NGC 5273 


) 


13.96 


0.65 


7.03 


1.32 


0.37 


14.50 


a 


NGC 5548 


n 


16.2! 


2.81 


0.65 


2.54 


0.20 


11.95 




NGC 5838 


\ 


14.84 


4.35 


2.57 


1.87 


0.93 




b,d 


NGC 5982 


n 


!2.69 


0.48 


1.73 


1.28 


0.06 


>21.0 




NGC 6340 


\ 


!2.54 


0.28 


2.46 


1.28 


0.59 


>20.0 




NGC 7457 


) 


13.46 


0.33 


2.32 


1.03 


0.35 


14.74 


a 


NGC 7626 


) 


13.24 


0.59 


1.84 


1.30 


0.36 


>19.0 


a 


NGC 7743 


\ 


11.59 


0.16 


5.36 


1.38 


0.50 


14.81 





Note. — Col. (1) Galaxy name. Col. (2) Profile type, where n — core, ) = intermediate, and \ = power law. Col. (3) 
Surface brightness at the break radius. Col. (4) Break radius. Col. (5) Transition parameter between outer and inner 
slope. Col. (6) Outer slope. Col. (7) Inner slope. Col. (8) -ff-band magnitude of the unresolved nucleus. Col. (9) Notes: 
(a) Lies in "gap region" 0.3 < 7 < 0.5. (b) Parameters for Nuker-law fit are highly uncertain due to extensive presence of 
dust. TTowever, the fit provides an estimate of the galaxy background for obtaining the nuclear magnitude, (c) Parameters 
for Nuker-law fit somewhat uncertain due to presence of strong edge-on disk, (d) Point-source magnitude of nucleus could 
not be determined. 
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TABLE 3: Comparison of the Nuker Parameters 



NICMQS WFPC1/WFPC2 



Galaxy jUf, rt q /3 7 jUt rf, q /3 7 Ref. 



NGC 221 


9 


.49 





.33 


4 


.66 


1 


.26 


0, 


.50 


12, 


,91 





.47 


1 


.39 


1 


,47 





.46 


1 


NGC 524 


13 


.73 


1 


.41 





.68 


1 


.69 


0, 


,03 


16, 


,12 





.32 


1 


.29 


1 


,00 





.00 


2 


NGC 3115 


13 


.38 


5. 


.39 


1 


.13 


1 


.80 


0, 


,73 


16, 


,25 


2 


.91 


1 


.47 


1 


.43 





.78 


2 


NGC 3379 


12 


.80 


1 


.58 


1 


.82 


1 


.45 


0, 


,18 


16, 


,14 


1 


.74 


1 


.59 


1 


.43 





.18 


2 


NGC 3900 


12 


.35 





.23 





.29 


1 


.66 


0, 


,51 


16, 


,28 





.56 





.12 


1 


.61 





.60 


3 


NGC 4278 


12 


.80 





.97 


1 


.63 


1 


.39 


0, 


,02 


15, 


,98 





.89 


1 


.45 


1 


,31 





.00 


4 


NGC 4406 


12 


.98 


1 


.00 


3 


.31 


1 


.16 


0, 


,00 


16, 


,08 





.94 


4 


.13 


1 


.04 





.04 


4 


NGC 4472 


13 


.52 


2 


.63 


1 


.89 


1 


.29 


0, 


,04 


16, 


,66 


2 


.41 


2 


.08 


1 


.17 





.04 


2 


NGC 4589 


12 


.05 





.21 


1 


.09 


1 


.18 


0, 




16, 


,61 





.75 





.43 


1 


.62 





.00 


4 


NGC 4636 


14 


.60 


3 


.44 


1 


.69 


1 


.56 


0, 


,13 


17, 


,73 


3 


.21 


1 


.64 


1 


.33 





.13 


2 


NGC 5982 


12 


.69 





.48 


1 


.73 


1 


.28 


0, 


,06 


15, 


,58 





.47 


2 


.15 


1 


,18 







4 


NGC 6340 


12 


.54 





.28 


2 


.46 


1 


.28 


0, 


,59 


15, 


,98 





.54 


1 


.73 


1 


.24 





.71 


3 


NGC 7626 


13 


.24 





.59 


1 


.84 


1 


.30 


0, 


,36 


16, 


,25 





.40 


1 


.53 


1 


.22 





.00 


4 



References.— (1) Lauer et al. 1998; (2) Faber et al. 1997; (3) Carollo & Stiavelli 1998; (4) Carollo et al. 1997a. 
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(this work) 



Fig. 2. — Comparison of aperture photometry derived from NICMOS and ground-based images. The measurements have been made with 
~10" diameter apertures. The solid line denotes equality. The average difference between the two sets of data is (Am^) = {■mjj{HST) — 
m/f(ground)) = — 0.006±0.11 mag. 
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Fig. 3. — Dependence of the inner cusp slope 7 on (a) absolute magnitude of the galaxy (M^ ) and (6) breaJj radius (r;,). The distribution 

of 7 is shown in panel (c). Filled symbols represent core galaxies, and open symbols denote power-law galaxies. The galaxies from this study 
are plotted as circles, those from Faber et al. (1997) are shown as triangles, and those from Rest et al. (2001) as squares. Objects in our 
sample which fall in the "gap region" (0.3 < 7 < 0.5) appear as asterisks, while those from Rest et al. (2001) are shown by plus signs. Two 
galaxies which fall in this region from Faber et al. (1997) are shown by crosses. 
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Fig. 4. — Comparison of the values of (a) 7, (6) /3, and (c) r;, derived from NICMOS images to those derived from optical images. The solid 
line denotes equality. Filled circles correspond to NGC 524, NGC 4589 and NGC 7626, which are known to have dust in the central regions 
as seen on optical images. The central profiles of these objects most likely were depressed in the optical by dust extinction, leading to 7(opt) 

Kl 0. 
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Fig. 5. — Central-parameter relations for core galaxies. The core parameters (r'f, and fif,) are well correlated with the central stellar velocity 
dispersion ((tq) and, with greater scatter, with the total -B-band absolute magnitude {M^ ). Panels (a) and (b) plot only data from this 
study (circles), whereas panels (c) and (d) include measurements from Faber et al. (1997; triangles) and Rest et al. (2001; squares). The 
outliers NGC 404 and NGC 4636 are labeled. 
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Fig. 6. — Fundamental-plane relation between ^n, and for core galaxies. The outliers NGC 404 and NGC 4636 are labeled. 
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Fig. 7. — Correlation of nuclear point-source magnitude with (a) extinction-corrected Ha luminosity and (6) central stellar velocity 
dispersion. Filled symbols denote nuclei in core galaxies, open symbols represent nuclei in power-law galaxies, and asterisks mark objects 
with intermediate inner slopes (0.3 < 7 < 0.5). Upper limits are indicated with arrows. Two objects with uncertain nuclear profiles appear 
as triangles. 
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